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I.  INTRODUCTION 


The  majority  of  tactical  missiles  used  against  surface 
targets  are  SEA-SKIMNING  and  fly  as  close  as  possible  to  the 
surface  of  the  water  in  order  to  make  difficult  the 
detection  and  reaction  of  the  enemy. 

This  work  addresses  the  problem  of  analysing  the  design 
of  the  altitude  control  system  of  a  supersonic  SEA-SKIMMING 
missile  in  both  aspects  of  SENSITIVITY  and  ROBUSTNESS,  with 
respect  to  the  variations  of  the  aerodynamics  parameters. 

Those  variations  can  develop  from  uncertainties  in  the 
model  as  well  as  from  changes  in  the  flight  conditions  of 
the  vehicle. 

The  SENSITIVITY  analysis  makes  it  possible  to  verify, 
for  a  given  model  which  parameters  are  more  important  for 
the  desired  flight  path. 

The  ROBUSTNESS  analysis  has  the  objective  of  verifying 
if  the  system  could  be  affected  by  perturbations  or  noise 
and  will  be  the  base  on  calculations  of  modified  feedback 
gains  which  are  capable  of  producing  the  desired  flight  of 
the  missile  despite  perturbations. 

The  contents  of  the  different  sections  are  as  fallowing: 

Section  II  presents  the  model  and  the  state  equations  of 
the  system  as  well  as  the  calculations  to  complete  the  chain 
of  control  of  the  missile,  considering  the  restrictions  of 


the  SEA-SKIMMING  scenario,  where  the  altitude  control  system 
cannot  tolerate  any  overshoot  to  avoid  that  the  missile  hits 
the  water. 

Section  III  develops  the  SENSITIVITY  ANALYSIS  taking 
into  consideration  the  variation  o-f  the  eigenvalues  and  time 

response  o-f  the  system  o-f  state  equations  for  variations  o-f 

the  aerodynamic  parameters  in  a  range  o-f  +25X  o-f  the 

original  value. 

In  section  IV  is  described  the  ROBUSTNESS  ANALYSIS  based 
upon  the  minimum  singular  value  o-f  the  return  di-f-ference 
matrix  as  function  of  the  frequency;  section  V  introduces  a 
design  technique  using  optimization  routine  in  order  to 
obtain  a  "complete  robust"  solution  and  section  VI  presents 
the  final  conclusions. 


A.  THE  MISSILE 

The  missile  will  be  considered  as  presented  by  Arrow 
CRef.lI.  This  model  has  been  selected  in  order  to  avoid 
classification  problems.  The  geometry  as  well  as  its 
characteristics  are  presented  in  Appendix  A. 

The  states  that  were  considered  in  our  model  are  defined 
in  Table  I  as  well  as  the  relevant  aerodynamics  parrmsters. 

In  the  Figure  I I. 2  we  have  the  schematic  representation 
of  the  uncoupled  pitch  control  chain  of  the  missile, 
neglecting  the  yaw  and  roll  movements; 


Figure  II. 2  -  Block  Diagram  of  the  Pitch  Control  Chain  of 
the  Missile. 


where 

T)  .  -  "bias"  acceleration; 

zcb  ’ 

-  commanded  acceleration; 

-  achieved  acceleration; 
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v(t)  -  velocity; 

H(t)  -  altitude; 

H  —  commanded  altitude, 

s 

The  Trans-fer  Function  G<s)  represents  the  missile  itsel-f 
and  the  other  blocks  are  the  controller.  The  input  to  this 
block  will  be  the  "commanded"  acceleration  and  the  output 
will  be  the  achieved  acceleration. 

The  commanded  acceleration  comes  from  the  feedback  of 
the  states  of  the  controller  (compensated  state,  velocity 
and  altitude)  added  to  a  "bias”  acceleration. 

In  the  diagram  of  Figure  II. 3,  we  have  a  complete 

representation  of  the  system  with  the  states  used  in 
modelling  our  missile. 

The  "STATE  EQUATIONS"  which  represent  our  system  are  the 
f ol lowing: 


><^=  -150.0  -  2,646.66  x  705.75  x^ 

x,=  1.353  +  X,  -  1.353  ij  +  1.334  H 

2  13  'zc  s 

><3=  -6.572  1  “  5  X  3  ^».572  -6.481 

><„=  -44.332  X  -  59.109  x 

H  S  fc 

X _=  X  -  0.1482  X  -  0.0395  x  (II.  1) 

S  H  5 

x,=  -188.4  X  +  188.4  x. 

o  »  S 
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-0.4608  x^-  2.231  x^  -  0.3406  x^  +  2.231  x^  + 

-15.0949  X,  -  20.13  x^  -  0.143  x,  +  0.4608  i) 
s  6  7  'zc 

-  0.4544  H 


X  .=  -50.0  X 

S  8 


^9  ^10 


495.0  X ,  +  495.0  H 

9  5 


x^„=  -17.644  X,  -  4.705  x 
10  s  tf 


This  system  o-f  equations  can  be  represented  in  matrix 
-form,  where  we  have: 


X  «  A  X  +  B  u 


-150  0 


O  2,646.66  -705.74  0000 


1.35  0  1.0  O 


-6.5  O  -5.0  O 


0  0  0  0 


0  0  0  0 


0  -44.33  -59.11  0  0 


0  1.0  -0.15  -0.04  0  O 


-188.4  188.4  0 


-0.46  -2.2  -3.4  2.2  -15.1  -20.13  -0.14  0 


O  0 


O  0 


0  O 


0  O 


O  -50.  -495  0 


0  0  0  1 


0  -17.64  -4.71  0000 


F  = 


B 


O  O 

-1.353  1.334 

6.572  -6.481 

O  O 

O  O 


O  O 

O  O 

0.461  -0.454 

O  495.0 


O  O 

u  =  -  F  X 


0  0  0  0  0  O  0 

O  0  0  0  O  O  O 


0.049  0.98615  1,404 
0  0  0 


The  vector  input  "u"  has  two  components 
-  the  commanded  acceleration;  and 


H  -  the  desired  altitude, 
s 

The  input  is  considered  as  two  parts;  one 

corresponds  to  the  feedback  of  the  states  (from  the 
controller)  and  another  is  the  "bias"  acceleration  of  Ig. 
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TABLE  I 

DESCRIPTION  OF  THE  SYSTEM  STATES  AND 
AERODYNAMIC  COEFFICIENTS 


X 

X 

X 


1 


2 


3 


X 

X 

X 

X 

X 

X 

X 


H 

s 

6 

7 

8 

9 

t  0 


.  -filtered  acceleration 

.  acceleration  compensator 

.  acceleration  compensator 

.  q  -  pitch  angular  rate 

. . .  <x-  angle  of  attack 

. .  &-  pitch  tail  incidence 

P 

commanded  pitch  tail  incidence 

.  x^-  controller  compensator 

. . .  H<t)  -  altitude 

.  v(t)  -  velocity 


AERODYNAMIC  COEFFICIENTS 


C 

C 

C 

C 


slope  of  curve 

mo: 


,  change  in  C  per 
mSp  ’  m 


slope  of  curve  of 
No: 


Np  change  in  C^^  per 


of  pitching  moment  coefficient 

degree  pitch  control  incidence 

normal  force  coefficient  C,,x  oc 

N 

degree  pitch  control  incidence 
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The  model  can  be  separated  in  two  main  parts;  the 


CONTROLLER  and  the  AUTOPILOT/AIRFRAME. 

B.  THE  CONTROLLER 

The  controller  was  developed  based  in  the  procedure 
introduced  by  Dowdle  CRe-f.2D. 

The  missile  was  considered  as  a  point  mass  and, 
therefore  the  altitude  as  -function  o-f  time  is  calculated  by 
double  integration  o-f  the  achieved  acceleration. 

The  controller  was  developed  without  take  in 
consideration  the  missile  itsel-f,  i.e.,  the  trans-fer 
function  G(s)  in  Figure  II. 2  was  assumed  as  "unity". 

According  to  the  [;Ref.2],  this  procedure  is  valid  if  the 
bandwidth  of  S(S)  is  larger  than  the  crossover  frequency  of 
the  controller. 

Due  to  the  constraints  of  the  scenario,  the  controller 
has  to  be  designed  in  such  way  that  no  overshoot  will  be 
tolerate.  Further  requirements  are  that  the  angle-of-attack 
and  the  acceleration,  have  to  remain  below  a  certain  level. 

The  controller  are  responsible  for  the  commanded 
acceleration  by  feedback  of  its  states. 

The  technique  used  to  solve  the  feedback  problem  was  the 
LINEAR  QUADRATIC  STATE  FEEDBACK  REGULATOR.  This  solution  has 
a  large  stability  margin  with  respect  to  phase  variations, 
as  well  infinite  gain  margin. 


2 


The  "cost  -function"  to  be  minimized  is  the  -folloMing: 

a> 

J  =  r  C  zTt)  a  z(t)  +  ult)  R  u(t>  }  dt  (II-2) 
0 

In  Figure  1 1. 4  we  have  the  schematic  representation  o-f 
the  CONTROLLER. 


Figure  1 1. 4  -  The  Controller 


The  LEAD  compensator  indicated  in  the  controller  is  used 
in  order  to  improve  the  transient  response  characteristics. 
Its  trans-fer  -function  is  : 


G 

c 


(s) 


K*(s-»-w  ) 
z 

(s+w  ) 
p 


where  K  =  10.0 

w  =  0.5  rad/sec 
z 

w  =  50.0  rad/sec 
P 

In  order  to  have  the  state  equations  o-f  the  CONTROLLER, 
we  draw  in  FIGURE  II. 5  its  signal  -flow  graph  using  the 
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adequate  transformation  of  the  compensator  and  showing  the 
states  that  will  be  feedback. 


H 

s 

f 


Figure  II. 5  Signal  Flow  Graph  of  the  Controller 

The  state  vector  is  defined  as: 

x<t)  =  C  x(t)  (H(t)-H  )  v(t)  3  ^  (II. 2) 

c  s 

In  the  Figure  1 1.5  the  states  of  the  controller  are 
represented  as  indicated  in  Table  I. 

From  the  diagram,  we  have  the  following  STATE  EQUATIONS: 

x=Ax+Bu  (II. 3) 

z  =  D  X 
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The  state  matrices  are: 


D  =  C  1  K  O  D 

Since  we  consider  to  feedback  all  the  states,  let  in 
equation  1 1 . 1 

Q  =  I  (Identity) 

and  R  =  0.04 

Using  the  OPTSYS  program  to  solve  that  OPTIMUM  CONTROL 
problem,  the  following  feedback  gains  were  obtained; 

=  0.049 

=  0.98615 
F^  =  1.4044 

correspondents  to  the  feedback  of  the  states  indicated 
in  Figure  II. 5. 
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C.  THE  AUTOPILOT/AIRFRAME 

The  controller  was  developed  assuming  the  transfer 
function  of  the  Autopilot/Airframe  as  unity. This  assumption 
implies  that  the  achieved  acceleration  is  equal  to  the 
commanded  acceleration;  however,  they  are  different  due  to 
the  aerodynamic  of  the  missile. In  order  to  take  that  in 
consideration  we  have  to  introduce  the  autopi 1 ot /ai rf rame 
model . 

Normally,  the  sea-skimming  is  a  SKID-TO-TURN  missile 
(turns  instantaneously) ,  but  due  to  problems  of 
classification  we  have  used  for  our  Sea — Skimming  model  the 
Autopilot  given  by  Arrows  CRef.ll  for  a  BANK-TO — TURN 
missile  (has  to  bank  in  order  to  turn)  that  was  designed  to 
fly  at  30,000.0  ft  instead  of  sea  level. 

Take  in  consideration  that  the  main  objective  of  this 
work  is  to  establish  a  procedure  for  sensitivity  and 
robustness  analysis  and  that  we  are  studying  only  the 
vertical  movement  of  the  missile  uncoupled  of  any  other 
movement, we  consider  that  the  approximation  is  useful. 

In  Figure  II. 6  we  have  the  block  diagram  of  the  missile, 
representing  the  plant  without  the  controller. 

The  state  equations  derived  from  this  diagram  are  on  the 
Appendix  A. 
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qfrad/sec) 


D.  THE  TIME  RESPONSE 

Using  the  software  CONIROLS,  the  system  represented  by 
equations  <11. 1)  was  simulated  and  the  time  response  for  the 
altitude,  velocity,  acceleration  and  angle  of  attack  were 
obtained. 

The  altitude  is  in  Figure  II. 7  and  we  can  see  that  the 
missile  arrives  at  the  desired  altitude  in  about  10  seconds. 

No  overshoot  is  shown  and  a  small  steady-state  error  is 
presented  characterizing  our  Type  O  system. The  steady-state 
error  was  less  than  5’/.. 

The  angle-of-attack  , reproduced  in  Figure  I I. 8,  varies 
between  3.7  deg  and  -1.2  deg  satifying  the  condition 
required  by  the  airframe,  as  indicated  in  CRefs.  1  and  23. 

The  acceleration  is  in  Figure  II. 9  and  the  velocity  is 
in  Figure  II. 10  representing  a  smooth  behavior  of  the 
system. 

In  the  next  chapter  will  be  presented  the  Sensitivity 
Analysis  using  the  developed  model. 
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Figure  11.8  -  Ang 1 e-of -Attack (deg )  vs.  Time(sec) 
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Figure  II.  10  -  Veloci  ty  <-f  t /sec )  vs.  Ti(ne(sec) 
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III.  SENSITIVITY  ANALYSIS 


A.  INTRODUCTION 

The  objective  o-f  the  SENSITIVITY  ANALYSIS  is  to 
analyse  how  the  behavior  o-f  the  system  is  af-fected  by 
variations  of  the  aerodynamic  parameters. Those  variations 
are  due  to  uncertainty  o-f  the  design  and/or  problems  ocurred 
during  the  flight. 

From  the  developed  model  we  can  see  that  the  important 

parameters  to  our  system  are  C  ,  ,  C  and  C^,  ,  as 

^  ^  moc’  Noc’  m&  N&’ 

described  on  TABLE  1. 

The  variations  could  be  interpreted  as  perturbations  in 
our  STATE  EQUATION  SYSTEM  (A-matrix)  and  the  analysis, 
basically,  take  into  consideration  the  change  in  the  poles 
positions  due  to  change  in  the  parameters. 

As  we  have  feedback,  we  replace  the  A  matrix  by  the 
augmented  A-matrix  <  A^  =  A  -  B  F  ) . 

The  sensitivity  of  a  matrix  where  the  elements  are 
time-invari ant ,can  be  verified  in  different  ways.  Among 
those  procedures,  some  are  enumerated  below. 

First,  using  the  method  introduced  by  Golub  and  Van  Loan 
CrBf.33  we  can  analyse  the  matrix  to  verify  how  it  is  "ill- 

conditioned";  where  an  "ill-conditioned"  matrix  means  a  very 
sensit’ve  system. 
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In  order  to  accomplish  the  analysis,  we  have  to 
calculate  the  "condition  number".  It  is  de-fined  as: 

K(A^)  =  II  I  I  II  aI^ I  I 

o  o  a 

where 

I  I  A^ I  I  is  the  norm  of  augmented  A  and 

II  A  ^11  is  the  norm  of  inverse  augmented  A 

di 

The  condition  mumber  depends  upon  the  considered  norm 
but  it  can  be  proven  that  the  condition  of  the  matrix  will 
be  the  same  for  all  norms. 

When  we  use  the  2-norm  one  has 

K  (A  )  =  *  ^^a^  (III. 1) 

a  (A  ) 

-  a 

where  s  (A  )  and  «  (A  )  are  the  maximum  and  minimum 
a  -  a 

singular  value  of  the  A-augmented  matrix,  respectively. 

Calculating  the  condition  number  for  the  different 
values  of  the  parameters  it  is  possible  to  establish  how  the 
sensitivity  of  the  system  changes  when  there  are  changes  in 
some  elements  of  the  state  matrix. 

The  problem  of  this  kind  of  analysis  is  that  we  can  not 
verify  the  sensitivity  of  the  system  with  respect  to  a 
specified  element  of  the  state  matrix  but  only  the 
sensitivity  of  the  overall  system. 

In  Frank  CRef.4],  we  have  another  method  where  the 
analysis  is  done  based  upon  the  STATE  SENSITIVITY  EQUATIONS 
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that  could  be  characterized  in  time  domain,  frequency  domain 


orin  terms  of  a  performance  index. 

The  basic  idea  is  to  define  "sensitivity  functions"  and 
analyse  the  behavior  of  the  system  by  those  functions. 

The  analysis  starts  with  the  system  defined  by: 

X  =  f  (X,t,U,o£)  (III. 2) 

where: 

-  ^  represents  the  state  vector; 

-  t  is  time  ; 

-  u  is  the  input  vector  ;  and 

-  o{  represents  a  vactor  of  the  different  parameters. 

If  ^  is  changed  by  Aot  (some  perturbation)  and  using 
Taylor's  expansion, we  have: 

Ax  (t,oc)  ^x  j  A  <x.  (III. 3) 

aoc.  ^ 

The  subscript  q£  indicates  that  the  partial  derivative 

0 

is  taken  at  the  nominal  parameter  values. 

The  partial  derivative  is  called  the  trajectory  vector  ( 

X  )  with  respect  to  the  j—  parameter  and  has  the  same 


dimension  as  the  state  vector. 

The  components  of  that  vector  are  defined  as  the 
TRAJECTORY  SENSITIVITY  FUNCTION. 


X  ,  .  (  t  ,oC„  ) 
1  J  ’-0 


dx  (t  ,QC) 
da. 


(III. 4) 
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The  equation  III. 4  represents  the  partial  derivative  of 

the  i—  state  with  respect  to  the  j—  parameter.  Assuming  we 
have  n  states  and  m  parameters,  it  will  be  possible  to  form 
a  n  X  m  matrix  called  TRAJECTORY  SENSITIVITY  MATRIX  (III.5). 
The  matrix  is  also  called  the  JACOBIAN  MATRIX. 


(III. 5) 


The  change  in  the  behavior  of  the  system  will  be 
determined  by: 


r 

^  =  X  (t,o(„)  =)  X.  Aot .  (III. 6) 

’-0  -  Z_.  J  J 

j  =  i 

The  analysis  is  based  in  how  the  sensitivity  function 
changes  when  we  change  the  parameters.  The  result  could  be 
visualized  by  plotes  of  X  as  function  of  time,  as  states  of 
the  expanded  system. 

Finally,  we  consider  the  EIGENVALUE  SENSITIVITY 
presented  by  Porter  &  Crossley  CRef.5]. 

This  latter  procedure  was  choosen  to  be  used  in  the 
analysis  of  this  thesis,  where  the  sensitivity  of  the  sistem 
is  related  to  the  variation  of  the  eigenvalues  of  the 
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augmented  matrix  (A^  =  A  -  B  F)  with  respect  to  the 

variation  of  the  parameters. 

In  order  to  further  develop  our  analysis,  the  concept  of 
MODAL  ANALYSIS  will  be  introduced. 

B.  THE  MODAL  ANALYSIS 

The  concept  is  based  in  generating  the  input  vector  of  a 
system  by  linear  feedback  of  the  state  vector  in  such  way 
that  the  prescribed  eigenvalues  are  associated  with  the 
dynamic  model  of  the  resulting  close-loop  system. 

As  an  example,  let  us  assume  a  scalar  system: 

><  =  a  X  (t)  +  b  u(t)  ;  (III.  7) 

if  u<t)  =  0  ====>  X  =s  X  „  exp  (at) 

where  exp  (at)  defines  the  "mode"  rf  the  system. 

If  a  ^  0  the  system  is  stable  and  if  a  ^  0  the  system  is 
unstable. 

Assuming  g  as  the  feedback  gain,  we  have: 

u(t)  =  g  X (t)  (111.8) 

Inserting  the  value  of  u  from  equation  1 1 1. 8  in  equation 
III. 7: 

X  =  (a  +  b  g)  X  (t)  (III. 9) 

and,  solving  for  x(t); 

x(t)  =  X (0)  expC(a  +  bg)tl 

Now  the  mode  of  the  system  is  defined  by  the  exp(a+bg). 
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When  we  have  a  state  vector,  the  system  is  represented 
by 

x=Ax+Bu  (III. 10) 

Considering 

u  =  —  F  X 

and  assuming  all  the  states  as  observable  and 
controllable,  the  substituting  value  o-f  u  in  equation 
( 1 1 1 . 10) ,  we  have 

x=(A-BF)x  or  x=A  x  (III. 11) 

a 

The  precise  nature  o-f  the  -free  motion  o-f  the  continuous 

time  system  -following  any  disturbance  can  be  described  in 

terms  o-f  eigenvalues  and  eigenvectors  of  the  augmented  plant 

matrix  (A  ) . 

a 

Assuming  A  has  n  distinct  eigenvalues  (X,,X_,...,X  ); 

then  it  also  has  n  corresponding  linearly  independent 

eigenvectors  u  -u  ,...,u  such  that 
18  n 

A  u.  =  X.  u. ;  i=l,2,...,n  (III. 12) 

ai  11^  ’ 

Using  the  calculated  eigenvalues  and  eigenvectors  we  can 

write  the  modal  matrix  of  A  as: 

a 

U=Cu  u  ...ul 
18  n 

where  each  column  of  U  is  the  eigenvectors  of  A^. 

Based  in  the  modal  matrix  U,  the  equation  (III. 12)  could 
be  written  as 
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A  U 
a 


A  U 


( 1 1 1 . 13) 


with  A  as  a  diagonal  matrix  o-f  the  eigenvalues. 

The  matrix  A  is  obtained  from  equation  (III. 13): 

U“‘  A  U  =  A  (III. 14) 

a 

In  addition  to  the  eigenproperties  of  A  ,  the 

a 

corresponding  properties  of  its  transposed  matrix  (A')  play 

d 

an  important  role  in  the  modal  analysis. 

The  matrix  A'  has  the  same  eigenvalues  of  A  ,  but  a 

different  sets  of  eigenvectors.lt  can  be  represented  as 

A'  V.  =  X.  V.:  j  =  1,2,...  ,n 

a  j  J  J*  ’  ’  ’ 

where  v^  represents  the  eigenvectors  of  A^. 

In  CRef.51  it  is  shown  that 

V  u  .  =  u  .'  V  .  =  .  ( 1 1 1 .  16) 

J  1  1  J  iJ 

where  .  is  the  Kronecker  delta;  i.e. 

1 J  ’ 


^iJ 

=  1 

if 

II 

and 

j 

=  0 

if 

i  ^  j 

Thus  in  matrix  form 

V'U  =  I  (III. 17) 

where  the  collumns  of  V  are  the  vectors  v. . 

1 

Using  the  modal  matrix  of  A^  and  introducing  a  new  state 
vector  y(t) ,  we  have: 

X (t)  =  U  y (t)  (III. 18) 
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The  derivative  o-f  both  sides,  gives  us: 


X  '  t )  =  U  y  ( t ) 


(III. 19) 


wi  th 


X  (t)  =  A  X  (t> 
a 


(II  I. 20) 


Combining  equations  ( 1 1 1 . 18) , ( 1 1 1 .  19)  and  (III. 20),  we 


have: 


y  =  U  ^  A  U  y(t) 


and  using  the  equation  (III. 14), 


y(t)  =  A  y(t) 


(III. 21) 


Because  the  matrix  A  is  diagonal  the  equations  are 
uncoupled  and  the  solutions  o-f  the  system  represented  by 
(III. 16)  are  given  by 


y. (t)  =  y. (0)  exp(x.t) 

i  1  1 


(1 1 1.22) 


Aplying  equation  III. 18, 


(t)  =  Cu^  u^  ...  u^:r  y^ 


(I  II. 23) 


which  implies  that 

x(t)  =  u,y,  (0)  Bxp(A,t)  u_y_(0)  exp(x_t)  + 


Taking  the  value  o-f  y(t)  in  equation  (III.  18), 


y (t)  =  U  ^x (t) 


(III.24) 


Using  equation  (III. 17), 


y .  (0)  =  V X  (o) 
'^1  1 


(III. 25) 
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Combining  equations  (III. 23)  and  (III. 25): 

x(t)  x(0)  U.V.'  exp(X.t>  (III.  26) 

j/  111 

Equation  I I I. 26  shows  that  the  system  represented  by 
(II 1.8)  has  its  modes  described  by  the  eigenvectors  of 

and  its  transpose  and  eigenvalues  of  matrix.  Also,  in 

order  to  have  the  system  stable,  x (t)  - >  0  as  t  - >  ®  , 

which  happens  if  and  only  if 

Re  <  0  ,  for  all  the  values  of  i 

In  the  Modal  Control  CRef.51  all  this  procedure  is  shown 
for  the  case  when  some  of  the  eigenvalues  are  equal  which 
does  not  apply  to  our  system. 

As  the  role  of  the  eigenvalues  is  so  important,  we  can 
verify  how  the  sensitivity  of  our  control  1  system  by 
analysing  the  change  in  the  eigenvalues  of  the  matrix  A  or 
the  augmented  A  for  the  close  loop  system  when  any  element 
is  changed. 

In  order  to  avoid  calculating  the  new  eigenvalues  for 
each  different  A,  the  equation  (1 11.26)  can  be  used  to  find 
approximated  values  that  can  be  analysed  analytically. 
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C.  EIGENVALUE  SENSITIVITY 

In  our  system  Me  consider  the  change  in  the  elements  of 
the  matrix  due  to  variations  in  the  aerodynamic 
coefficients. 

Basically,  Me  concerned  Mith  the  eigenvalues  that  are 
changed  in  the  direction  of  the  right  hand  plane, making  the 
system  unstable  or  affecting  the  time  response. 

The  result  of  the  analysis  could  be  plots  or  tables 
relating  the  eigenvalues  Mith  the  variation  of  the 
parameters. 

For  each  change  in  the  parameters.  Me  have  a  neM  matrix. 

Approximations  to  the  eigenvalues  associated  Mith  the 
different  plant  matrices  in  the  vicinity  of  the  original  one 
could  be  used  Mithout  having  to  compute  the  eigenvalues  of 
each  neM  matrix. 

We  describe  the  approach  presented  in  the  reference 
CRef.5]  applying  the  "first  order  eigenvalue  sensitivity"  to 
the  change  of  one  of  the  parameters  and  doing  the  complete 
analysis  Mith  the  calculated  eigenvalues. 

1 .  First-order  Eigenvalue  Sensitivity 
From  equation  III. 11,  Me  have: 

X  (t)  =  A  X  (t) 
a 

Let 's  assume 


Mhere  a^^^  is  any  element  of  A^. 
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Differentiating  equation  (111.12)  with  respect  to 


a, ,  indicates  that 
kl 


u^  +  A  3ui  =  ^Xi  +  X  ^ui 

da,  ,  da,  ,  da,  ,  da,  , 

kl  kl  kl  kl 


<111.27 


multiplying  (III. 27)  by  v'  then  gives 

v.'  dA  u.  +  V.'  A  dui  =  v*  dxi  u.  +  X.  v'  dui 

1  -r-  1  1  T -  1  T -  1  11  T — 

da.  ,  da.  ,  da.  ,  da 


kl  — kl  ““kl 

As  A  is  equal  to  x^^  ,  we  have 


kl 


V.'  dA  u.  =  v'  dxi  u. 
kl  kl 

Which  reduces  to  the  set  of  equations 

=  V.  (k)  u.(l)  ;  i=k=l  =  l  ,2, .  .  .  ,n  (111.28) 

^"kl  ' 

Due  to 


1^  =  ^k  ^jl 

^®kl 


In  equation  (111.28),  v^(k)  and  u^(l)  represent  the 

k—  element  of  v^  and  the  1—  element  of  u^  ,  respectively. 

Those  coefficients  may  be  considered  as  the 
elements  of  a  set  of  n  eigenvalue  sensitivity  matrices. 

The  new  eigenvalues  (x^)  of  the  altered  matrix, 

when  the  element  a^^  is  changed  will  be  calculated  as 

X .  =  X .  +  _ _i  A  a,  , 

i  1  T -  kl 
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Mhen  mb  have  more  than  one  parameter  altered,  the 
linearity  could  be  used  and  the  -final  eigenvalue  will 
obtained  -from  the  summation  o-f  the  sensitivity  coe-f f icients 
multiplied  by  the  respective  change  in  the  coefficients. 


As  an  example  of  this  procedure,  let's  consider  our 

system,  changing  the  value  of  C  such  that 

msp 

AC,  =0.02 

msp 

When  changed  as  specified,  the  following 

elements  of  the  augmented-A  matrix  of  the  state  equations 
presented  in  the  Chapter  II  are  changed: 


A  (4,6):  Aa„^  =  14.778 

A  H 

A  (7,6) :  Aa,^  =  5.035 

The  eigenvalues  and  eigenvectors  of  the  augmented  A 
matrix  and  of  its  tranpose  were  calculated  using  the  program 
Mat lab. 


a)  Eigenvalue  : 


X^  =  -159.75  +  j  18.85 

A  X  =Cv  (4)  u  (6)  1  A  a  +  Cv  (7)  u  (6)  1  A  a 

Xa  X  vO  X  X  rO 

v^(4)  =  -0.0148  -  j  0.0039 
v^(7)  =  1 

u^(6)  =  0.014  -  j  0.027 
X^  =  -159.815  j  18.85 
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b)  Eigenvalue  ; 

^2  =  -159.815  -  j  18.85  (complex  conjugate  o-f  x^) 

c)  Eigenvalue  X^  ; 

Xj  =  -50.0 

AX  =  Cv  (4)  u  (6) :  Aa^^  +  Zv,i7)  u,(6)]  Aa 

^  i  a  Ho  3  3 

Vj(4)  =  0.0003 
Vj(7)  =  0.0039 

Uj(A)  =  -0.001 

ys 

Xj  =  -49.99998 

ys 

d)  Eigenvalue  x^  : 

=  -8.61  +  j  7.91 

AX^  =  i:v^(4)  u^(6)3  Aa^^  +  :v^(7)  u^  (6)  3  Aa^^ 
v^(4)  =  -0.1667  -  j  0.1322 
v^(7)  =  1 

u^(6)  =  -0.0334  -  j  0.132 
X^  =  -8.95  +  j  7.6355 

/N 

e)  Eigenvalue  x^  : 

Xj  =  -8.95  -  j  7.6355 

/S 

f)  Eigenvalue  x^  : 
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=  -1.82  +  j  2-02 

AX^  =  u^(6)D  (6)  3  Aa 

v^(4)  =  0.1718  -  j  0.0161 
v^(7)  =  -0.5737  +  j  0.1689 
11^(6)  =0.01  +  j  0.0363 

=  -1.846  +  j  2.05 

g)  Eigenvalue  x^  : 

X^  =  -1.846  -  j  2-05 

/S 

h)  Eigenvalue  x^  : 

X^  =  -2.71 

AX  =  Cv  (4)  u,(6>3  Aa^^  +  Cv„(7)  u  (6)3  Aa 
v^ (4)  =  0. 1741 
v^<7)  =  -0.6597 
Uj(6)  =  -0.0396 

X,  =  -2.681 

i>  Eigenvalue  X^  : 

X^  =  -0.5 

AX^  =  CVy(4)  u^(6)3  Aa^^  +  Cv^(7)  u^(6)3  Aa 
V ,(4)  =  0. 1718  +  j  0.016 

5 


v^(7)  =  -0.5737  -  j  0.1689 

u^<6)  =  0.01  -  j  0.0363 

X.  =  -0.4741 

9 

j)  Eigenvalue  x  ^  : 

X^o  =-0.14 

(4)=  0. 1741 

1  0 

V  (7)=  -0.6597 
1  0 

u, „ l6)=  -0.0396 
1  0 

X, .  *  -0.1104 

1  0 

For  the  analysed  parameter,  our  system  is 
unsens’ ti ve,  taking  in  consideration  the  small  dif-ferences 
between  the  new  and  original  eigenvalues. 

Repeating  this  procedure  -for  all  dif-Ferent  parameter 
values  will  give  us  the  approximation  of  the  new  eigenvalues 
and  that  can  be  helpfull  in  a  qualitative  analysis  of  what 
parameter  has  more  effect  over  our  model . 

As  we  want  to  verify  the  effect  in  the  time 
response  and  due  to  the  Software  available  the  analysis  will 
be  done  completely  with  the  calculated  exact  values. 


46 


Using  the  Controls  Program,  we  will  verify  the  variation 
of  the  eigenvalues  of  the  augmented  A  matrix  (  A  -  FB)  with 
respect  to  the  variation  of  the  parameters  and  by  simulation 
we  analyse  the  effects  on  the  time  response. 

At  the  end  of  the  analysis  we  compare  the  results  of  the 
example  used  for  the  first  order  approximation  with  the 
exact  values  that  were  calculated. 

As  reference,  the  eigenvalues  of  the  original  system 

are: 

=  -159.746  +  j  18.979 

Xj  =  -159.746  -  j  18.979 

=  -49.9992 

X^  =  -8.60697  +  j  7.90968 

Xj  =  -8.60697  -  j  7.90968 

x^  =  -1.81854  +  j  2.0238 

x^  =  -1.81854  -  j  2.0238 

X,  =  -2,70266 

X.  =  -0.503819 

9 

X  „=  -0.143073 
1  0 


Each  aerodynamic  coe-f -f  icient ,  that  takes  part  o-f  the 
considered  model  Mas  changed  up  to  +2Sy.  and  the  eigenvalues 

of  the  perturbed  system  matrix  ( A-augmented )  were  analysed 
as  well  as  the  TIME  RESPONSE  of  the  system  . 
a  •  O 

mcx 

□n  TABLE  II  we  have  the  eigenvalues  related  to  the 

considered  values  of  C  and  as  we  can  see,  with  exception 

mo( 

of  the  7—  eigenvalue,  the  changes  are  small  and  the  system 
is  pratically  insensitive  to  that  variation. 

The  time  responses  are  shown  in  Figures  III.l,  111.2,111.3 
and  I I I . 4  . 

-  Angl e~of -attack  (Figure  III.l)  -  changes  in  the  same 

direction  of  C  with  the  minimum  value  changing  about 
mo«  ^ 

0.5  degrees  in  the  negative  direction.  The  maximum 

value  stays  almost  constant  ; 

-  Acceleration  (Figure  III. 2)  -  the  change  was  less  than 

2 

5  ft/sec  either  the  maximum  or  the  minimum  values; 

-  Velocity  (Figure  III. 3)  -  almost  the  same  ;  and 

-  Altitude  (Figure  III. 4)  —  almost  the  same. 


The  change  in  the  eigenvalues  are  on  TABLE  III  ;  the 

time  response  is  on  Figure  111.5,111.6,111.7  and  III. 8. 

The  variation  of  the  eigenvalues  are  slightly  higher 

with  changes  in  C  but  the  system  stays  stable  and  the  time 

m&p 

response  has  a  small  variation. 
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-  Angle-o-f-Attack  (Figure  1 1 1. 5)  -  The  minimum  values 

increases  negativelly  about  0.5  degrees  but  in  the 
contrary  direction  o-f  the  coefficient,  i.e.,  increasing 
the  value  of  C  negativelly,  the  minimum  was  smaller: 

-  Acceleration  (Figure  III. 6)  -  small  changes; 

-  Velocity  (Figure  I I I. 7)  -  small  changes; 

-  Altitude  (Figure  III. 8)  -  pratically  the  same. 


c. 


C 


No; 


The  change  in  the  eigenvalues  are  on  TABLE  IV  and  the 
time  response  is  on  Figures  111.9,111.10,111.11  and  III. 12. 

The  changes  in  the  eigenvalues  are  small  but  the  effect 
on  the  time  response  is  more  sensible,  mainly  on  the  angle- 
of -attack . 

-  Angl e-of -Attack  (Figure  I I I. 9)  -  the  maximum  and  the 

minimum  increase  1.0  degree  in  the  positive  and 
negative  directions,  respecti vel 1 y ; 

-  Acceleration  (Figure  III. 10)  -  maximum  changes  of  lOX 
in  the  negative  direction  (minimum)  and  20'/.  in  the 
positive  direction  (maximum) . These  values  are  smaller 

when  we  increase  the  value  of  C.,  ; 

Noc 

-  Velocity  (Figure  III. 11)  -  small  changes; 

-  Altitude  (Figure  III. 12)  -  does  not  change. 


The  change  in  the  eigenvalues  are  on  Table  V  and  the 
time  response  are  on  the  Figures  111.13,111.14,111.15  and 
III.  16 

-  Angl e-of -Attack  (Figure  III. 13)  -  Changes  almost  1.0 

degree  with  the  maximum  and  minimum  increasing  when  the 
value  of  the  parameter  increases; 
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Acceleration  (Figure  111.14)  -  the  maximum  increases 

with  the  value  o-f  the  parameter; 

Velocity  (Figure  111.15)  -  pratically  the  same  -for  all 
the  considered  values  of  the  parameter; 


Altitude  (Figure  III. 16) 


constant . 


TABLE  II 


C  -  SENSITIVITY 
moc 


-0. 

04 

-0.  05 

-159.741 

J 

18.9843 

-159.743 

+ 

j 

18.9816 

-159.741 

-  j 

18.9843 

-159.743 

- 

j 

18.9816 

-49.9992 

-49.9992 

-8.55822 

+  j 

6. 37 130 

j 

7. 18736 

-8. 55822 

-  j 

6.37130 

-8.56386 

- 

j 

7. 18736 

-2.31740 

+  j 

2.95574 

-2. 14117 

+ 

j 

2.43275 

-2. 31740 

-  j 

2.95574 

-2. 14117 

- 

j 

2.43275 

-1.79944 

-2. 14200 

-0.51650 

-0.510282 

-0. 14304 

-0. 14306 

-0. 

07 

-0.  08 

-159. 748 

+  j 

18.9762 

-159.751 

+ 

j 

18.9735 

-159. 748 

-  j 

18.9762 

-159.751 

- 

j 

18.9735 

-49.9992 

-49.9992 

-8.66197 

+  j 

8.55389 

-8.72178 

+ 

j 

9. 14021 

-8.66197 

+  j 

8.55389 

-8.72178 

+ 

j 

9. 14021 

-1.47981 

+  j 

1.81088 

-1.22077 

+ 

j 

1 . 70081 

-1.47981 

-  j 

1.81088 

-1.22077 

- 

j 

1 . 70081 

-3. 27131 

-3.67125 

-0. 49779 

-0. 49148 

-0. 14308 

-0. 14308 

5 


TABLE  III 


C  -  SENSITIVITY 

m&p 


c  ^ 

m  &p 

-0. 

06 

-0. 

07 

-162.604 

+  J 

16.7131 

-161.200 

+ 

j 

17.9141 

-162.604 

-  j 

16.7131 

-161.200 

- 

j 

17.9141 

-49.9992 

-49.9992 

—5. 74603 

+  j 

8. 86120 

-7. 13033 

j 

8.50992 

-5. 74603 

-  j 

8.86120 

-7. 13033 

- 

J 

8.50992 

-1.24694 

+  j 

1.66681 

-1.54230 

+ 

j 

1 . 80335 

-1.24694 

-  j 

1.66681 

-1.54230 

- 

j 

1.80335 

-3.86370 

-3.30548 

-0.  49103 

-0.49828 

-0. 14311 

-0. 14309 

C  , 
m&p 

-0. 

09 

-0. 

10 

-158.243 

+  j 

19.9192 

-156.678 

+ 

j 

20.7492 

-158. 243 

-  j 

19.9192 

-156.678 

- 

j 

20. 7492 

-49. 9991 

-49. 9990 

-10. 1621 

+  j 

6.98960 

-11.7920 

+ 

j 

5.53519 

-10. 1621 

-  j 

6.98960 

-11.7920 

- 

j 

5.53519 

-1 . 47981 

+  j 

1.81088 

"2.03752 

j 

2.55411 

-1.47981 

-  j 

1.81088 

-2.03752 

- 

j 

2.55411 

-2.27254 

-2.02276 

-0.50786 

-0.51120 

-0. 14307 

-0. 14306 
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TABLE  IV 


C.,  -  SENSITIVITY 

Noc 


'Not 


'No£ 


0. 

13 

0. 

14 

-159.750 

+ 

j 

18.9895 

-159. 748 

+ 

j 

18.9842 

-159. 750 

- 

J 

18.9895 

-159.748 

- 

j 

18. 9842 

-49.9992 

-49.9992 

-8.76667 

+ 

j 

8.27012 

-8.68784 

+ 

j 

8.09494 

-a. 76667 

- 

8.27012 

-8. 68784 

- 

j 

8.09494 

-1. 18629 

+ 

j 

1.74321 

-1.47163 

+ 

j 

1 . 84903 

-1. 18629 

- 

j 

1.74321 

-1.47163 

- 

j 

1 . 84903 

-3.63456 

-3.22775 

-0.48878 

-0.49663 

-0. 14314 

-0. 14310 

0. 

16 

0. 

17 

-159.744 

+ 

j 

18.9736 

-159.741 

+ 

j 

18.9683 

-159. 744 

- 

j 

18.9736 

-159.741 

- 

j 

18.9683 

-49.9992 

-49.9992 

-8.52441 

+ 

j 

7.71274 

-8.44077 

+ 

j 

7.50213 

-8.52441 

- 

j 

7.71274 

-8. 44077 

- 

j 

7.50213 

-2. 15317 

+ 

j 

2.32249 

-2.40105 

+ 

j 

2.67918 

-2. 15317 

- 

j 

2.32249 

-2.40105 

- 

j 

2.67918 

-2.20610 

-1.88565 

-0.51041 

-0.51651 

-0. 14305 

-0. 14303 
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TABLE  V 


Q,  -  SENSITIVITY 

N&p 


i  C..  0.02  0.03 

I  N&p _ 

-159.004 
-159.004 
-49.9996 
-9.37545 
-9.37545 
-2.24432 
-2.24432 
-1.98723 
-0.51356 
-0. 14304 


&p 

0. 

05 

0. 

06 

-160. 106 

+ 

j 

21.5632 

-160.459 

+ 

j 

23.8437 

-160. 106 

- 

j 

21.5632 

-160.459 

- 

j 

23.8437 

-49. 9990 

-49. 9988 

-8.25204 

+ 

j 

8.25703 

-7.91449 

+ 

j 

8.57438 

-8.25204 

- 

j 

8.25703 

-7.91449 

- 

j 

8.57438 

-1.57370 

+ 

j 

1.85345 

-1.35009 

j 

1.75196 

-1.57370 

- 

j 

1 . 85345 

-1.35009 

- 

j 

1.75196 

-3. 18783 

-3.60970 

-0.49848 

-0.49275 

-0. 14310 

-0. 14312 

12.1138 

-159.748 

+ 

j 

18.9842 

12. 1138 

-159.748 

- 

j 

18.9842 

-49.9992 

7. 1 1519 

-8.68784 

+ 

J 

8.09494 

7. 11519 

-8.68784 

- 

j 

8.09494 

2.56109 

-1.47163 

+ 

j 

1 . 84903 

2.56109 

-1.47163 

- 

j 

1 . 84903 

-3.22775 

-0.49663 

-0. 14310 

54 


D.  CONCLUSION 


As  a  conclusion  for  the  sensitivity  analysis,  we  have 
the  following: 

With  respect  to  the  eigenvalues,  the  pitch  control 

system  is  more  sensitive  to  changes  in  C  and  C  .  than 

msp  N&p 


changes  in  C  and  C 
mot 

Considering  the 
angle-of-attack  is 
higher  variation 


No,  * 

time  response, 

sensitive  to 

with  C.,  ;  the 
Not 


we  can  verify  that  the 
all  the  parameters  with 
acceleration  is  more 


sensitive  to  C  and  C^,  ;  the  velocity  has  the  maximum 

mot  Noc  ' 

value  changed  less  than  lOX  with  respect  to  all  the 
parameters;  and  the  altitude  is  insensible  for  the 
considered  variations. 

It  is  not  necessary  to  emphasize  the  importance  of  this 
analysis  for  the  design  of  the  control  system,  but  thinking 
only  in  the  eigenvalues  and  comparing  the  results  of  the 
first  order  approximation  presented  for  with  those  of 

Table  III  for  C  changed  to  -0.06  (6C  =  0.02)  we  can  have 

m  &p  m  &p 

an  idea  of  the  importance  of  the  parameter,  using  that 
approx i mat i on . 


55 


SENSITIVITY 


Figure  III.l  Angle-of  Attack  vs  Time 

■for  different  values  of  C 
_  mot 
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Figure  III. 3  Velocity  vs  Time 

for  different  values  of  C 

moc 
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Figure  III. 7  Velocity  vs  Time 

for  different  values  of  C 

msp 
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Figure  III. 9  Angle-o-f-Attack  vs  Time 

•for  di-fferent  values  o-f  C.. 
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Figure  III. 13 


Angl e-of -Attack  vs  Time 

for  different  values  of  C. , , 
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Figure  III. 14  Acceleration  vs  Time 

■for  different  values  of 


69 


TIME(SEC) 


-  SENSITIVITY 


TIME(SEC) 


SENSITIVITY 


TIME(SEC) 


IV.  ROBUSTNESS  ANALYSIS 


A.  INTRODUCTION 

In  the  previous  chapter  we  take  into  consideration  the 
change  in  the  eigenvalues  and  in  the  time  response  o-f  the 
system  due  to  variations  on  the  aerodynamic  parameters  that 
take  part  into  the  model;  under  the  robustness  analysis  we 

1 

! 

veri-fy  what  we  can  expect  o-f  the  system  when  it  is  affected  j 

by  any  kind  of  perturbation. 

i 

I 

The  main  idea  is  to  use  the  robustness  analysis  to  j 

design  a  new  set  o-f  -feedback  gains  in  order  to  have  the  same 
behavior  in  di-fferent  situations  and/or  environments  which 
act  as  perturbations  to  our  missile. 

In  classical  -frequency  domain,  -for  single-input  single- 
output  systems (SISO) ,  a  robust  design  can  be  achieved  using 
Bode,  Nyquist  or  Nichols  plots.  With  these  techniques  it  is 
possible  to  de-fine  gains  that  give  us  gain  and  phase  margins 
■for  a  "robust"  system.  For  multi-input  multi-output  systems 
(MIMO) ,  the  classical  techniques  are  no  longer  valid. 

Taking  this  in  consideration,  we  have  analysed  the 
robustness  o-f  the  original  design  by  means  o-f  the  minimum 
singular  value  of  the  return  difference  matrix  for  different 
frequencies  and  we  will  use  the  same  technique  to  improve 
the  control  system. 
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In  order  to  show  why  we  are  applying  this  technique,  we 
start  with  the  Nyquist  Criterion  -for  a  SISO  system  and 
extend  the  analysis  to  a  MIMO  system,  using  the  minimum 
singular  value  theory,  as  explained  by  Lehtomaki ,  Sandell 
and  Athans  in  ERef.61. 

B.  THE  NYQUIST  CRITERION 

Let's  take  in  consideration  the  SISO  system  represented 
in  Figure  IV. 1 . 


Figure  IV. 1  Single-Input  Single-Output  System 

where  G(s)  includes  the  controller  and  plant  dynamics 
with  unity  -Feedback  . 

The  Nyquist  criterion  states  that  i-f  the  open  loop 
trans-fer  -function  G(s)  does  not  have  any  pole  in  the  right 
half  s  plane,  then  the  locus  of  G(s)  will  not  encircle  the 
point  (-1,0)  in  the  Nyquist  plot,  where  jw  is  substituted 
for  s  and  the  axes  represent  the  real  (Re)  and  imaginary 
(Im)  parts  of  G(Jw)  for  various  frequencies. 
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A  Nyquist  Plot  is  given  in  Figure  IV. 2  and  it  can  be 


used  to  design  feedback  gains  that  ensure  a  robust  system. 


The  measure  of  the  system  stability  are  Gain  and  Phase 
margin  obtained  from  the  Figure,  as  indicated. 


Figure  IV. 2  Nyquist  plot  -  Stable  System 

The  Phase  Margin  (PM)  is  the  additional  phase  lag  at  the 
gain  crossover  frequency  required  to  bring  the  system  to  the 
limit  of  instability  and  the  Gain  Margin  (GM)  is  the 
reciprocal  of  the  magnitude  |G(jw)|  at  the  frequency  where 


the  phase  angle  is  -ISO  .  The  gain  crossover  -frequency  is 
the  frequency  at  which  |G(jw) I  is  unity. 

The  Phase  and  Gain  margins  are  a  measure  of  how  close  a 
polar  plot  is  to  the  -1  jO  point. 

Considering  the  system  subject  to  additive  perturbation, 
as  we  have  in  the  diagram  of  Figure  IV. 3. 


Figure  IV. 3  Additive  Perturbation 

In  order  to  have  the  system  stable,  despite  the 
perturbation,  we  need 

|A6(jw)  I  <  II  +  G(jw)  I  (IV. 1) 

As  we  can  see  in  Figure  IV. 2,  this  condition  ensures  a 
stable  system. 

This  idea  can  be  extended  to  MIMD  problems  through  the 
use  of  matrix  norms  and  of  applying  the  Multivariable 
Nyquist  Theorem. 
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c. 


MULT I -INPUT  MULT I -OUTPUT  SYSTEM 


The  generalization  o-f  the  SISO  theory  discussed  in  the 
previous  section  has  been  made  -for  the  MIMO  problem. 

Let's  consider  our  system  as  represented  in  chapter  III 
and  apply  the  transformation  in  order  to  have  the  open  loop 
transfer  function  between  input  and  output. 

X  =  A  X  +  B  u  (IV. 2) 

Applying  Laplace  Transform  to  the  equation  (IV. 2),  we 

have 

sX  =  AX  +  BU  (IV. 3) 


where  X  and  U  represent  the  Laplace 

u(t),  respectively. Taking  the  value  of  X 

X  =  Csl  -  Al"^B  U 
As  Y  =  C  X 

We  have  Y  =  C  Csl  -  Al'^B 

U 


Transform  of  and 

in  equation  (IV. 1) 

(IV. 4) 

(IV. 5) 

( IV. 4) 


This  transfer  function  corresponds  to  G(s)  in  Figure 

IV. 4. 


Figure  IV. 4  Return  Difference 
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From  Figure  IV. 4  it  is  possible  to  de-fine  the  "Return 
Dif-ference  Matrix"  (Kwakernaak  and  Sivan,  CRef.63)  that 
corresponds  to  opening  the  loop  at  a  given  point  and 
connecting  an  external  input  variable  at  this  point;  -for 
example,  at  point  1.  Assuming  the  input  u  as  zero,  we  have: 

Y^(s)  =  -  G(s)  H(s)  V  (s)  (IV. 5) 

where  V  (s)  corresponds  to  the  new  input  variable  and 
Y^(s)  the  "returned  variable". 

The  difference  between  Yj(s)  and  V(s>  is 

V(s)  -  Y^(s)  =  C  I  +  G(s>H(s)l  V(s)  (IV. 6> 

The  matrix  I  G(s)H(s)  is  defined  as  the  "Returned 

Difference  Matrix",  with  G(s)  equal  to 

G(s)  =  CCsI  -  A]"S  (IV. 7) 

As  explained  by  Lehtomaki ,5andel 1  and  Athans  in  CRef.7], 
the  Multivariable  Nyquist  Theorem  is  derived  from  the 
rel ati onshi p 

detCI  +  G(s)H(S)D  =  $cl (s)  (IV. 9) 

lol  (s) 

where  ^1  (s)  =  det(sl  -  A)  corresponds  to  the 

characteristic  polynomial  of  the  open  loop  transfer  function 
and  lol (s)  =  det(sl  -  A  +  BF)  corresponds  to  the 

characteristic  polynomial  of  the  closed  loop  transfer 
function  from  the  system  represented  in  Figure  IV. 4. 


The  multivariable  Nyquist  theorem  which  requires  that  a 
closed  loop  stable  system  have  the  same  number  o-F 
counterclockwise  encirclements  of  the  origin  by  the  locus  of 
the  determinant  of  I  +  G(s)H(s)  as  the  number  of  open  loop 
poles  that  are  unstable. 

If  I  +  G(s)H<s)  is  quasi -si ngul ar ,  a  small  change  in  G 
may  make  the  matrix  singular  ;  this  causes  the  detCI  + 
G(s)H(s}]  to  become  zero  and  the  Nyquist  encirclement  count 
to  change  indicating  an  unstable  system. 

Basically, in  order  to  analyse  the  robustness  of  the 
system  we  have  to  verify  how  close  the  return  difference 
matrix  is  to  being  singular  as  a  function  of  frequency  (u»). 

The  natural  measurement  of  the  singularity  is  the 
minimum  "singular  value",  since  this  is  the  tightest  norm. 

The  singular  value  of  a  matrix  M  is  defined  as 

5.  =  €X.  (  (IV.  10) 

1  1 

M  represents  the  conjugate  transpose  matrix  of  a 

|_| 

generic  M  matrix  and  any  eigenvalue  of  the  product  of  f1 
times  M. 

If  the  minimum  singular  value  is  close  to  zero  the 
matrix  is  quasi -si ngul ar  and  the  system  is  not  robust. 

If  one  assumes 

»<I  +  G(jw>H(jw)>  >  (IV. 11) 
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then  it  can  be  shown  CRe-f.  71  that  the  gain  and  phase 
margins  -for  the  system  may  be  represented  by  equations 
(IV. 12)  and  (IV. 13),  i.  e. , 

Qain  margin:  GM  =  1  (IV. 12) 

Phase  margin:  PH  =  +  cos  ^C1  -  (IV. 13) 

2 

Universal  gain  and  phase  margins  curves,  presented  by 
Gordon  CRe-f.  81  and  repeated  in  Figure  IV. 5  allow  the 
designer  to  pick  a  singular  value  that  corresponds  to  a 
desired  stability  margin  for  the  system  to  be  designed. 

From  singular  value  plots  the  we  will  be  able  to 
identify  the  critical  frequencies  or  range  of  frequencies 
where  the  singular  value  is  below  the  chosen  level  from  the 
Figure  IV. 5. 
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D.  SINGULAR  VALUE  ANALYSIS 

Using  the  program  developed  by  Gordon  CRe-f.Sl,  with 
small  changes  to  make  it  compatible  with  the  size  o-F  our 
system  and  to  input  the  matrices  more  easily,  we  calculate 
the  minimum  singular  value  o-f  the  return  difference  matrix, 
with  the  frequency  varying  from  0  to  50  rad/sec. 

The  program  utilizes  the  IMSL  subroutines  to  generate 
the  return  difference  matrix  as  well  as  to  made  the  singular 
value  decomposition,  given  the  minimum  singular  value  for 
different  frequencies. 

The  return  difference  matrix  depends  on  the  point  where 
we  open  the  loop  represented  in  Figure  IV. 5.  i-or  the 
explained  case,  considering  point  1,  we  have  the  "output" 
difference  matrix,  that  corresponds  to  I  +  G(s)H(s>;  if  the 
loop  IS  open  at  point  2  we  have  the  "input"  case  and  the 
return  difference  matrix  will  be  I  +  H(s)G(s). 

From  the  Figure  IV. 5  we  choose  as  reference,  for  our 
analysis,  the  minimum  singular  value  as  0.6  that  gives  us  a 

0 

phase  margin  of  about  35  to  characterize  a  "robust"  system. 

The  minimum  input  singular  value  is  plotted  in  Figure 
IV.  6. 

It  shows  that  the  system  is  robust,  except  for 
frequencies  between  0.6  and  5.0  rad/sec. 
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The  results  related  to  the  minimum  output  singular  value 
are  in  Figure  IV.  7  and  the  robustness  o-f  the  system  is  poor 
•from  0  to  18  rad/sec. 

Those  results  indicate  that  the  behavior  of  the  system 
under  perturbations  should  present  problem  for  frequencies 
below  10  rad/sec. 

To  confirm  this  expectation,  we  have  plotted  the  BODE 
diagram  (magnitude)  for  the  closed  loop  system  considering 
the  angl e-of -attack  and  altitude  with  respect  to  both  inputs 
(acceleration  and  desired  altitude). 

The  diagrams  are  on  Figures  IV. 8  to  IV. 11,  where  can  see 
the  attenuation  at  the  low  frequencies  indicating  that  we 
can  not  expect  good  response  at  those  frequencies  and  one 
could  further  expect  problems  with  robustness  as  is  evident 
from  the  singular  value  plots. 

On  the  next  chapter  we  try  to  .mprove  the  system 
designing  a  new  set  of  feedback  gains,  in  order  to  have  a 
robust  system. 
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Figure  IV. 9  Bode  Magnitude  -  Angl e-of -Attack  vs. 
Desired  Altitude 
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FREQUENCY  -  RADIANS  PER  SECOND 
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FREQUENCY  -  RADIANS  PER  SECOND 
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Figure  IV. 11  Bode  Magnitude  -  Altitude  vs. 
Desired  Altitude 
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FREQUENCY  -  RADIANS  PER  SECOND 


V.  IMPROVING  THE  DESIGN 


Under  this  chapter  we  try  to  improve  the  system  by 
calculating  new  -feedback  gains  that  will  yield  a  more  robust 
system,  taking  into  consideration  the  results  indicated  in 
Figures  IV. 6  and  IV. 7. 

The  POP’  AR  program,  developed  by  Gordon  CRe-f.Gl  will  be 
used.  A  numerical  optimization  technique  is  applied  in  order 
to  increase  the  minimum  singular  value  o-f  the  considered 
return  dif-ference  matrix,  there-fore,  a  resultant  robustness 
o-f  the  design. 

The  development  in  this  chapter  will  be  preceded  by 
comments  on  numerical  optimization  and  a  description  of  the 
computer  program. 

A.  OPTIMIZATION 

The  optimization  was  accomplished  by  means  of  the 
Automated  Design  Synthesis  Program  (ADS)  developed  by 
Vanderplaats  CRef.91. 

The  purpose  of  ADS  as  of  others  numerical  optimization 
routines  is  to  find  the  "best"  possible  solution  for  the 
problem,  starting  from  an  initial  set  of  variables  and 
updating  the  design  i teracti vel y .  The  problems  can  develop 
in  convergence  of  the  method  and  in  the  computer  time 
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needed.  I-f  the  problem  has  multiple  solutions,  the 
optimization  does  not  alMays  lead  to  the  absolute  optimum. 

The  ADS  program  is  designed  as  a  black  box  optimizer 
which  allows  the  user  to  choose  combinations  of  one 
dimensional  search,  optimization  algorithm  and  strategy. 

ADS  is  used  as  a  subroutine  and  the  parameters  that 
correspond  to  the  different  applications  are  chosen  by  the 
user,  as  explained  in  CRef.91. 

The  preferred  method  used  in  the  ADS  applications  is 
referred  to  as  Sequential  Unconstrained  Minimization 
Techniques  and  can  be  considered  as  a  method  that  starts 
with  an  objective  function  and  the  constrains  combined  into 
an  augmented  objective  function  and  then  minimizing  this 
function  as  if  the  problem  was  unconstrained. 

ADS  employs  penalty  function  techniques  as  well  as  an 
Augmented  Lagrange  Multiplier. 

As  presented  in  CRef.83,  the  iteration  between  the 
user's  program  and  the  ADS  routine  can  be  represented  by  the 
block  diagram  of  Figure  V.l;  where  the  "info"  parameter  is 
used  as  a  "flag"  for  a  dialogue  between  the  user  and  the  ADS 
routines. 

The  details  for  utilization  of  ADS  are  presented  by 
Vanderplats  in  CRef.Sl. 
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Figure  V. 1  Organization  of  ADS  Program 
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B.  POPLAR  PROGRAM 


The  Gordon's  program  has  the  capability  of  calculating 
the  feedback  gains  that  improve  the  response  of  the  system 
either  by  pole  placement  and/or  by  increasing  the  minimum 
singular  value  of  the  return  difference  matrix. 

The  main  objective  is  to  improve  the  robustness  by 
increasing  the  minimum  singular  value  above  a  desired  level 
chosen  from  the  universal  curve  reproduced  in  Figure  IV. 5. 

"Optimum"  values  of  the  feedback  are  calculated  by 
minimization  of  an  Objective  function  whose  pole  placement 
part  IS 


1 


OBJ 


i  =  l 


£ 


where 

X  -  real  part  of  the  desired  eigenvalue  i; 
di 

X  -  real  part  of  the  computed  eigenvalue; 

R  . 

I 

X  -  imaginary  part  of  the  desired  eigenvalue  i; 
■^di 

Xj  -  imaginary  part  of  the  computed  eigenvalue  i; 
i 


and,  for  the  minimum  singular  value  optimization,  we 

have 

OBJ  =2^  ^  ^  (jw,p)  3  J 

j 

where  indicates  the  desired  miriimum  singular  value 

and  ^  is  the  minimum  singular  value  at  a  certain  frequency. 
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The  optimization  procedure  changes  the  -feedback  gains 
until  the  minimum  singular  value  is  raised  to  the  desired 
1 evel . 

The  pole  placement  and  robustness  program  calculates  the 
return  di-f-ference  matrix  to  be  considered,  working  in  the 
complex  space  (as  the  analysis  is  conducted  in  the  frequency 
domain)  and  the  objective  functions,  calling  the  ADS  routine 
to  make  the  optimization. 

Input's  for  the  program  are  the  matrices  of  the  state 
representation  of  our  missile  (  A  ,  B  ,  C  and  F) , where  the 
observation  matrix  (C)  was  specified  as  Identity  because  we 
are  assuming  all  the  states  are  observable  as  well  as 
control 1 abl e. 

Other  inputs  are  the  desired  minimum  singular  value, 
desired  pole  locations  and  frequency  interval  to  be 
considered. 

In  our  application  of  POPLAR  the  initial  values  of  the 
feedback  gains  are  those  of  the  original  system  and  the 
parameters  for  application  of  the  ADS  program  are; 

Strategy  -  Augmented  Lagrange  multiplier; 

Optimizer  -  BFGS  variable  metric  method  for  unconstrained 

mi ni mi zati on; 

One-Dimensional  Search  -  Polynomial  interpolation. 

Two  situations  will  be  considered;  the  INPUT  -  minimum 
singular  value,  where  the  return  difference  matrix  i s  I  +  HG 
and  OUTPUT-mi ni mum  singular  value  with  I  +  GH  as  the  return 
difference  matrix. 
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C.  INPUT  -  MINIMUM  SINGULAR  VALUE 

The  program  was  -first  run  to  improve  the  minimum 
singular  value  -for  the  input  situation  within  the  same 
interval  o-f  -frequencies  used  in  Chapter  IV. 

Due  to  the  large  CPU  time  involved,  the  program  was  run 
-for  the  critical  band  o-f  -frequencies  and  the  result 
evaluated  for  the  complete  interval. 

The  minimum  singular  value  taken  as  reference  from  the 

0 

Universal  curve  was  0.6  that  implies  a  phase  margin  of  35  . 

We  have  an  objective  function  taking  into  consideration 
the  improvement  of  the  minimum  singular  value  as  well  as  a 
part  corresponding  to  the  pole  placement.  The  main  goal  is 
to  improve  the  singular  value  but  with  a  solution  such  that 
the  poles  are  in  positions  where  the  time  response  satisfies 
the  requirements. 

The  weight  of  the  pole  placement  part  was  considered  as 
10’/.  of  the  singular  value  part. 

The  program  was  started  considering  only  part  of  the 
feedback  matrix  as  free  parameters,  increasing  this  number 
until  a  reasonable  solution  was  found. 

The  best  solution  was  found  keeping  the  original 
feedback  gains  and  calculating  gains  to  feedback  all  tb-? 
other  states  to  generate  the  commanded  acceleration. 

The  computer  output  of  the  best  solution  with  the 
corresponding  inputs  is  given  in  Appendix  B. 
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The  -feedback  gains  that  improve  the  robustness  of  the 
system  are  the  following 

0.04  -1.15  0.03  -1.7  -9.07  0.32  7.07  0.049  0.986  1.4 
000  0000  O  0  0 

The  plot  of  the  minimum  singular  value  for  different 
frequencies  is  given  in  the  Figure  V.2,  along  with  the 
values  for  the  original  system. 

As  we  can  see  the  minimum  singular  value  was  increased 
in  the  critical  values  of  frequency  with  the  minimum 
changing  from  0.40  to  0.63  that  indicates  an  improvement  in 

0  0 

the  phase  margin  from  20  to  35  .  Also  the  bandwidth  where 
the  system  was  less  robust  was  changed  from  5.0  rad/sec  to 
approximately  1.0  rad/sec. 

The  robustness  of  the  control  system  for  all  the 
considered  frequencies  is  better,  arriving  close  to  0.9  with 
a  considerable  improvement  at  low  frequencies. 

Using  the  software  Controls,  the  time  response  was 
plotted  in  Figure  V.3. 

The  maximum  and  minimum  values  of  the  angl e-of -attack 
are  smaller  than  those  of  the  original  system;  the  maximum 
acceleration  and  velocity  are  slightly  lower  and  an 
undershoot  of  about  10'/.  appears  on  the  altitude,  that  could 
be  a  problem  in  some  flight  conditions. 
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We  now  turn  to  a  consideration  of  the  output  return 
difference  matrix  although  this  is  less  important  from  a 
practical  viewpoint. 
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In  the  Figure  V.4  we  have  plotted  the  output  minimum 
singular  values  that  correspond  to  the  run  where  the  input 
singular  values  were  improved,  as  well  as  the  original 
val ues. 

As  a  considerable  degradation  of  the  output  singular 
values  was  evident,  we  try,  in  the  last  part  of  this  work, 
to  arrive  to  a  situation  where  the  system  could  be 
reasonable  robust  in  both  cases,  increasing  the  output 
singular  value  even  with  some  reduction  at  the  input 
si tuation. 

After  several  runs  of  the  program,  using  different 
combinations  of  optimization  techniques  with  different 
starting  points,  it  was  verified  that  due  to  the 
characteristics  of  the  plant,  that  increases  in  the  input 
singular  values  resulted  in  decreases  in  the  output 
situation. 

In  order  to  obtain  some  improvement,  an  effort  was  made 
to  modify  the  system. 

The  first  step  was  to  look  at  the  controllability  matrix 
of  our  model ,  represented  by 


1.8556320-01 
2. 7256640+13 
2.7256640+13 
2.8611210+00 
2.2018110+01 


4.9481680+02 
-2.6904490+13 
-2.6904490+13 
-2. 8220090+00 
2. 1717100+01 


99 


-6.7817910+00 


6. 6891 110+00 


-2.2869750+00  2.2555630+00 

9.8554330+00  -9.7196140+00 

-1.5083820+01  1.4876850+01 

4 . 9054870-0 1  -4 . 2643980-0 1 

As  we  can  verify,  all  the  elements  are  non  zero,  this 
means  that  the  poles  can  placed  at  any  desired  position;  but 
a  close  look  at  the  2nd  and  3th  rows  indicates  very  high 
order  numbers  compared  to  the  others.  This  situation 
generates  some  numerical  difficulties  in  pole  placement. 

Using  engineering  judgement,  we  have  to  define  what  kind 
of  changes  have  to  be  made  for  achieving  a  more  robust 
system. 

First  we  try  to  get  a  balanced  A-matrix  by  changing 
units,  i.  e. ,  the  angles  will  be  in  radians  instead  of 
degrees  and  the  pitch  angular  rate  in  radians  per  second 
instead  of  degrees  per  second.  Despite  some  reduction  on  the 
numbers  was  not  possible  to  arrive  to  a  satisfactory 
solution . 

More  positive  effect  was  obtained  by  adding  the  effect 
of  the  angl e-of-attack  in  the  controller  (see  Figure  1 1.3). 

As  the  original  system  has  the  output  singular  values 
higher  than  0.5  for  frequencies  above  10  rad/sec,  we  have 
kept  the  original  feedback  gains  and  calculated  gains  to 
feedback  only  the  angl e-of-attack ,  pitch  rate  and  commanded 


actuator . 


The  gains  are  the  following: 


■  O  0  O  0.037  0.027  O  -0.02  0.049  0.986  1.404' 

F  = 

.000  0  00  O  0  o  o. 

The  result  is  presented  in  Figure  V.5,  compared  with  the 
original  singular  values  and  those  from  the  improvement  in 
the  input  case. 

The  output  singular  values,  compared  with  the  original 
was  slightly  improved  for  frequencies  below  5.0  rad/sec. 

Resultant  input  singular  values  are  plotted  in  Figure 
V.6,  and  they  are  lower  than  those  from  the  improvement  in 
the  input  situation  but  a  little  higher  than  the  original. 

The  time  response  were  practically  unchanged  with 
respect  to  the  original  model  and  it  is  presented  in  Figure 
V.7. 

Further  analysis  is  needed  to  obtain  a  solution  with 
high  output  singular  values. 
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OUTPUT  -  MINIMUM  SING.  VALUE 
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Figure  V.5  Output  Singular  Value 
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VI.  CONCLUSICNS 


The  sensitivity  analysis  based  in  the  change  in 
eigenvalues  represents  an  important  role  in  the  design. 

As  was  demonstrated  for  the  studied  missile  i'^s  behavior 
corresponds  to  a  system  insensitive  to  the  variation  of  the 
considered  aerodynamic  parameters. 

'^he  analysis  gives  a  good  indication  of  what  parameters 
have  to  be  precisely  determined. 

The  robustness  analysis  also  demonstrated  how  useful  the 
singular  value  analysis  is  as  an  auxiliary  tool  for  the 
designer . 

Singular  value  analysis  indicates  that  the  system  is 
robust  to  input  perturbations  but  is  deficient  in  robustness 
to  output  perturbations. 

The  physical  nature  of  the  problem  indicates  that  the 
prime  concern  should  be  with  respect  to  the  input 
perturbations;  the  output  perturbations  are  of  minor  concern 
but  should  be  kept  in  mind  if  unusual  conditions  should  be 
encountered  by  the  missile. 

□n  the  MIMO  design  a  commomly  used  method  is  the  Linear 
Quadratic  analysis  where  the  performance  levels  are  reached 
by  adjusting  weighting  terms  in  the  "cost"  function,  but  the 
results  of  that  method,  for  non-diagonal  R  matrices,  do  not 
necessarily  imply  a  robust  system.  Using  the  return 
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di-f-ference  matrix  we  can  improve  the  robustness  of  the 
design. 

One  advantage  of  this  procedure  is  that  it  permits  a 
high  level  of  interaction  betv^een  the  system  and  the 
designer. 

Further  problems  of  a  computer  nature  (i.e. ,  large  CPU 
time)  are  encountered  with  a  high  number  of  states.  Analysis 
must  be  confined  to  a  small  frequency  range  of  low  singular 
values  and  the  complete  system  return  matrix  singular  value 
are  calculated  with  the  determined  feedback  gains. 

One  aspect  that  should  be  considered  in  future 
development  is  to  improve  the  cpu  time  used  when  handle  a 
system  with  a  high  number  of  states. 
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APPENDIX  A 


MISSILE  DATA 

Under  this  Appendix  the  missile  sizing,  mass  properties 
and  aerodynamic  parameters  are  presented  as  given  by  Arrow 
CRe-f.l],  including  the  -figures. 

The  missile  is  1/6  scale  of  the  actual  circular  missile 
configuration  and  is  reproduced  in  Figure  A.l.  It  is  tail 
controlled  using  four  identical  control  surfaces  located 
0 

with  +30  dihedral. 

A.  BEQMETRY  AND  MASS  PROPERTIES 

In  table  VI  we  have  the  size  and  mass  properties,  with 
the  respactive  values,  used  for  development  of  the  state 
equations. 

Only  the  uncoupled  pitch  channel  was  considered, 
assuming  no  roll  movement. 


Table  VI 

MISSILE  GEOMETRY 

AND  MASS 

PROPERTIES 

Wei ght 

W 

2525 

1  bs 

Length 

1 

168 

i  n 

Di ameter 

d 

24 

in 

Reference 

Area  for  Coeff. 

S 

IT 

2 

ft 

Moment  of 

Inertia  about  y 

I 

804 

slug- 

-ft* 

yy 
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B.  AUTQPILOT/AIRFRAME 

Figure  A. 2  shows  the  block  diagram  of  the  autopilot  used 
to  develop  the  model  of  the  pitch  control  channel  of  the 
mi ssi 1 e. 

The  state  equations  that  correspond  to  the  autopilot  are 
the  following: 

X -150.0  +  150.0 

X  1 . 353  X  +  X  -  1 . 353  i) 

2  13  zc 

X  ,=  -6.572  X  -  5.0  X  +  6.572  i) 

3  13  zc 

x^=  -44.3316  Xg  -  59.11  x^ 

X  =  X  -  0.14«2  x^  +  0.0395 

S'?  5  o 


x=  -188.4  X  +  188.4  x_ 

w  O  / 

X  =  -0.4608  X  -  2.231  x  -  0.3406  x  +  2.231  x 

7  12  3  H 

-  15.095  X,  -  20.13  X,  -  0.1430  x  +  0.4608  t) 

S  6  7  ZC 


The  input  for  the  autopilot  is  the  commanded 
acceleration - ij 

zc 

The  values  related  to  the  aerodynamic  are  the  following: 


C 

mac 


C 


m&p 


C 


Noc 


O.  06 

0.08 

O.  15 

0.04 
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Dynamic  Pressure  -  q  -  1650 

Velocity  -  V  -  3025.46 


Ib/ft * 
■f  t/sec 


ill 
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q{rad/sec) 


APPENDIX  B 


The 

singular 


POPLAR  PROGRAM 

computer  output  where  the  improvement 
value  was  found  is  listed  below. 


in  the  input 


^iLCi  utrsfNo*  ovr  ai 


<«<<  HUH  •  soet  »>» 


WAX*  f.iooot  iNtT  pxeQ«  •.7ooo«  fxeo  srtf»  t.itoM 


HCtOMT  Im  l.tOO««MeiOHT  2«  •.•OCSttCfOKT  S« 


SVHINl*  l.tOOOl  9VniN0«  Ilf*  l.lftot  too*  • 

*••  thc  a  ^lakt  matxix 


•isp.oeoot 

9.00009 

•lOS.ltOOt 

9.99999 

9.90999 

«.««««# 

9.00009 

l.USU 

9.90009 

1.90009 

9.99999 

9.09999 

9.99009 

9.99999 

9.99999 

0.P00Q9 

9.09099 

•i.ijzu 

9.99099 

9.99099 

9.99000 

9.90099 

9.90999 

9.000P9 

9.90009 

9<OQOO« 

9.00009 

9.90099 

9.90909 

-44.11102 

•$9.11099 

9.09909 

9.99999 

9.Q0Q09 

9.00009 

9.90099 

9.90009 

9.90909 

1.99099 

•9.14929 

•9.919S9 

9.99999 

9.99909 

9.00099 

9.90009 

9.00009 

9.90009 

9.90099 

9.90009 

9.99009-199.19999 

I99.S9999 

9.99099 

0.00099 

9.00099 

*0.4<099 

•Z.2J099 

<.21099 

•IS. 99919 

-29.11009 

-9.14199 

9.90999 

9.90099 

9.90009 

9.00099 

9.99099 

9.99099 

9.99999 

9.99999 

9.99999 

9.99999 

-99.99099 

•099.90099 

9.90909 

9.00099 

9.90999 

9.99009 

9.99999 

9.99999 

9.99999 

9.99999 

9.90999 

9.00099 

1.90009 

9.00009 

9.00909 

9.99999 

9.99999 

•19.94499 

-4.79S99 

9.99999 

9.99999 

0.00009 

9.90009 

•••  TMf  • 

COHTROk  l^rtfT  MATRIX 

9.90099 

9.99909 

-I.ISSI9 

I.JS449 

*.sr:99 

-9.49199 

9.00009 

9.90999 

9.09099 

9.90909 

9.00009 

9.99999 

9.40099 

9.90009 

909. 99009 

9.90009 

9.90009 

9.00009 

9.90099 

•••  THC  C 

09UAVATI0H 

MATRIX 

• 

1.90099 

9.99909 

9.90999 

9.99999 

9.99999 

9.99999 

9.99999 

9.99999 

9.90099 

9.90099 

9.90099 

1.90009 

9.99999 

9.99999 

9.99999 

9.99999 

9.99999 

9.99999 

9.90999 

9.99009 

9.90999 

9.99009 

1.99999 

9.99999 

9.99099 

9.99999 

9.99999 

9.99999 

0.90009 

9.99099 

9.99999 

9.99999 

9.99999 

1.99999 

9.99999 

9.99999 

9.99999 

9.99999 

9.99099 

9.99009 

9.99999 

0.00990 

0.09099 

9.90909 

1.99909 

9.99099 

9.99909 

9.90009 
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FAOC  fft 


pii.ei  LATsmo*  OUT  ai 


0.00000 

0.00000 

e.oeoof 

0.00000 

o.oeeoo 

0.00000 

0.00000 

e.oeoeo 

1.00000 

0.00000 

0.00000 


THC 


0.00000 

0.90015 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 


0.00000 
1.40440 
0.00000 
0. 00000 


o.oeeoo 

0.00000 

0.00000 

o.eoeoo 

0.00000 


0.00000 

0.01000 


i.oooeo 
0.00000 
i. 00000 
0.00000 
o.oeoQO 


0. 00000 
0.00000 


0.00000 

0.00000 

0.00000 

o.oooeo 

0.00000 


t. 00000 

0.00000 

0.00000 

0.00000 

o.oeeoo 


0.00000 
t. 00000 
0.00000 
0.00000 
0.00000 


o.oooeo 

o.toooo 

l.ooeoo 

o.oooeo 

0.00000 


0.04900 

0.00000 


0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000 


P  FCeSOACK  MATRIX  ■■■ 


oRDCReo  corvLCx  ciocHUALucs  (i>^un 


139.00000 
139.00000 
•30.00000 
•0.00000 
•t.ooooo 
•9.00000 
•5.00000 
•I .00000 
•I .00000 
•I .00000 


19.09990 

19.09999 

0.00000 

T.ooooo 
0. 00000 
o.oooeo 
t.ooooo 
0.00000 
0.00000 


•••  INITtAC  PATA  15  ••• 


•  ••  T«  OROCRCO  COMPVrCO  eiOCNVACUCS 


-159.942]] 
•159.942)2 
•50.02094 
•0.41140 
-t. 43440 
•2.44024 
•1.04204 
-1.04204 
-0.50400 
-0.14S11 


-11.99300 

10.97300 

0.00000 

-7.09021 

7.09021 

0.00000 

-2.15090 

2.15090 

0.00000 

0.00000 


■  ■■  T>C  OROCRCO  ClOOIVtCTORO 


29.44)94 

-14,45100 

-29.44)94 

14.45100 

1.50)92 

0.00000 

-4.43)30 

11.39435 

-4.4)350 

-11.1945$ 

0.1)429 

o.toooo 

4.47929 

-1.12431 

4.47929 

1.32431 

5.00250 

0.90000 

1.79549 

o.foeoo 

0.24904 

0.095)0 

0.24904 

0.32102 

0.00000 

-1.94O00 

-|.)439l 

-1.04000 

1.14)90 

2.19104 

0.00000 

0.04739 

0.00130 

0.04739 

-I.OIIIt 

0.1322S 

0,00000 

0.00472 

0.00000 

-1.20790 

-0.44)41 

-1.20799 

0.44)4) 

-1.94491 

0.00900 

9.293T7 

0.44100 

1.29)77 

-4.44100 

-11.29992 

0.00000 

-4.42302 

-4.25209 

-4.42302 

4.25205 

-0.90271 

0.00000 

-01.01242 

0.00000 

-0.04045 

-•.)44)9 

-0.04045 

0.)4439 

•0.20442 

0.00900 

-7.73927 

-1.4:410 

—  9.7)927 

1.42410 

). 37471 

9.99000 

0.72477 

0.9)041 

0.72497 

•0.13041 

0.01012 

0.00000 

•0.00140 

9.90009 

0.00034 

0.00199 

0.00034 

-0.00197 

0.00351 

0.00000 

0.90900 

-0.34190 
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FruCi  LAT^ROt  OWT 


FAOC  90S 


A1 


e.9«ro0 

0.S4374 

-O.S04OS 

0.00000 

-0.11322 

0.002*9 

-0.21122 

-0.00299 

-0.2S444 

0.00000 

-0.24040 

0.00000 

-0.01491 

-1.00104 

-0.01491 

1.00104 

-0.2450S 

0.00000 

-1.14940 

-0.94011 

-l.SASii 

0.94010 

0.49479 

o.oeeoo 

t.tS044 

-0.01917 

0.21044 

0.01917 

0. 19194 

0.00000 

0.20211 

0.00000 

-O.lOSl? 

-0.19000 

-0.I0JI7 

0. 15000 

-0.17990 

0.00000 

-1.14140 

-l.S4S(0 

0.04451 

0.49030 

0.00000 

f.2SS04 

-0.02171 

0.21504 

0.02171 

0.1914S 

0.00000 

0.20190 

0.00000 

-0.04044 

-0.00944 

-0.04044 

0.00944 

-251.15575 

0.00000 

-1.20490 

-1.20494 

1 .4S904 

-12.59140 

0.00000 

2.72214 

•.040S4 

t. 72214 

-5.04014 

-140.74147 

0.00000  -1020.24701 

0.00000 

0.00004 

0.00010 

0.00004 

-o.eooii 

0.00042 

0.00000 

0.07750 

0.01127 

0.07750 

-0.OS327 

1.S2212 

f  .00000 

-i.24S9t 

-0.51411 

-0.24591 

0.51411 

14.00240 

0.00000 

104.95024 

f  .00000 

-O.OOS01 

-0.02004 

-0.00301 

0.02004 

-0.02110 

0.00000 

-0.40479 

-5.59004 

-0.40470 

0.09004 

-1.25274 

0.00000 

-0.44022 

1.40442 

-0.04022 

-1.40442 

-7.1)040 

0.00000 

-24.47534 

0.00000 

■  •a  T>C 

F  FCtOtACK 

NATRIK  ■■■ 

0.00000 

0.9R41S 

4.00400 

1.40440 

f .00000 

0.00000  0.00000 

0.00000 

0.00000 

0.04900 

0.00000 

0.00000 

0.00000 

0.00000 

0.01000 

0.00400  0.04000 

0.00400 

0.00400 

0.00000 

FRCOl^NCV 

HIM  AOO 

IN  SV 

2NO  AOO  IN  IV  MAX 

AOO  IN  SV 

MIN  RROO 

ADD  IN  SV 

0.70004 

0.55204 

1.25129  1.25329 

1.14110 

FRCQ 

SVAONO 

IVADXO 

svnm  ivNix 

SVTRtO 

SVMXO 

0.70004 

0.05114 

15.14094 

0.522202221.19727 

2.00400 

2.00400 

AAAAA  DODOOD 
AO  D 

AD  0 

AAAAAA  0  0 

AO  0 

AO  0 

A  A  OODDOO 


S1SSSS 

I 

S 

SSSS9 

s 

s 

ssssss 


RORTRAH  PROOAAH 
^00 

UTOMATKO  OISIOM  $VNTM»*l» 
V  C  K  S  I  0  N  1.00 


CONTROL  RARAMCrCRS 
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ritci  LATs^Noa  Out 


PAOC  #t4 


ISTRAT  ■  t  lOfT 

lOMAO  •  t  MOV 


I  ■  I  tfAlKT  •  !••• 

7  NCCN  •  $ 


OPTIMIZATION  RtSULTS 


OIJCCTIVC  FUNCTION  VA4.UV  t.iatT«C>t« 


DCSICN  VARUaLCS 


LOMca 
to  UNO 

-t.ttaeoc^os 
•o.iooooe*«« 
-9. iooooe>oi 
-9. )eoooe*9s 
•>9.19090e*9S 
-9, |0090e«99 
*9. 19090C«9$ 


VALUC 

9.i«>ire-9i 
•9.1t494e<91 
9.  J0T9Se-91 
•9.l««Sie*9t 
•  9.T9«MC*9t 
9.SI4t4e«99 
9.79l4Se«9t 


UFFOI 

tOUNO 

t.t9999e«9f 
9.  |999ae*95 
9.  i0900e«95 
9.iooeoe*9S 
9.1999ee*0S 
9.  i9oeoe*9i 
9.t9909e*99 


function  cvaluations  •  zii 
•  T»<  OFTlMlZet  CO»ruTCO  OVTFVT  11  ••• 


■  T»«  oADeiieo  cor^vrcD  eiocNVAcuei 


•  Mt.Tsni 
•ISt.TSIII 
•9C  .91:59 
•9.94499 
•9.94499 
-7.94445 
•5.99411 
•9.J74S4 
•9.17414 
•0.14105 


•19.94907 

19.94997 

9.90909 

7.97497 

9.99999 

9.90999 

0.17759 

9.00000 


•■t  THE  QFOCRCO  CIOCNVCCTORS 


27.97851 

•12.15019 

-27.97951 

12.15019 

1.71545 

9.99999 

•9.41152 

•9.74942 

-9.91152 

9.94992 

•2.91255 

9.99909 

•2.22249 

9.99999 

9.44219 

-9.57774 

9.49219 

9.57794 

-1.99572 

9.99999 

9.27174 

9.99291 

9.27174 

•9.09291 

9.14915 

9.99999 

9.19909 

•9.71511 

9.19999 

9.91511 

14.599  94 

9.99999 

•2.99941 

9.99099 

9.91279 

9.99427 

9.91299 

•9.99429 

•9.99719 

9.90999 

•1.11079 

•9.49259 

•1.13979 

9.49259 

•1.47471 

9.99909 

•9.71771 

9.19919 

-9.91771 

•9.19914 

•72.17257 

9.99099 

14.44447 

9.99999 

•9.19777 

9.91774 

-9.19979 

•9.91999 

•1.49999 

9.90999 

-9.94929 

•9.14999 

•9.99929 

9.14909 

•9.19921 

9.99999 

9.57172 

-9.47214 

9.50191 

9.49214 

•1.91192 

9.99999 

•9.99149 

9.99999 

9.92457 

9.99252 

9.92459 

-9.99252 

•9.99994 

9.99909 

9.99049 

9.99197 

9.99949 

•9.99197 

9.99577 

9.99909 

•9.99922 

9.94219 

-9.90022 

-9.99219 

9.27422 

9.90009 

9.17745 

9.99999 

*9.91291 

9. 94997 

•9.91291 

-9.94097 

9.24157 

9.90009 

-9.949SI 

-9.99219 

•9.44411 

9.74239 

•9.24514 

9.90099 

9.11517 

•9.94574 

9.11517 

9.04574 

•9.41442 

9.99009 

-9.29799 

9.99909 

9.92449 

•9.11947 

9.9:449 

9.41947 

•9.l959i 

9.90009 

-9. 10951 

•9.14151 

•9.19951 

9.14111 

-9.17497 

9.99999 

9.12744 

-9.94759 

9.12749 

9.94759 

-9.41911 

9.99999 

•9.29917 

9.99999 

9.92449 

•9.91044 

9.92449 

9.91944 

•9.19909 

9.99999 

-9.94995 

•9.44917 

•9.99995 

9.99437 

•241.94149 

9.99499 

9.17754 

-9.949C7 

9.19754 

4.94927 

9.99129 

9.99999 

•9.47299 

9.99999 

29.11471 

14.14297 

29.11491 

•14.14:97 

1792.99559 

9.99999 
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^ILCi  i^TSfRO*  OVT 


At 


fAOC 


t. 10095  i.oeoir  o.ooeos  -t.toet?  o.eoooo  •.•ooo«  ••.tooos  t.oo««« 

•f.eoofS  •«. 00490  -o.tsoot  i.ooeoi  -o.ootst  •.ooooo  -t.ozoor  -1.42104 

•2.02(47  1.421*4  •109.27l7t  t.lOOOt 

-9.00472  -0.027SS  •0.00472  •.027SI  -0.02205  0.00090  0.004SS  O.OKlt 

0. 0*435  -0. 01(15  0.2*554  0.00090  0.422(7  0.00090  0.1*520  1.517*0 

O.I*S2(  -1.SS7**  25.00540  0.09009 


T>«  0  PtCOOACK  HATRIX 


0.05*51  -1.14*44  0.05070  -1.(9552  -*.0(*74  O.SISS*  7.07442  0.04*90 

O.fOdS  1.40(40 

0.00900  0.00000  0.01000  0.00000  0.00000  0.00009  9.00009  0.90099 

9.00009  9.00009 


meauencv  mu  aoo  in  sv  2no  aoo  in  tv  tux  aoo  in  sv  nih  pnoo  add  in  sv 

0.70000  0.(1511  2. 0(2*5  2.0(295  1.1*025 

7XCQ  SVAOHO  SVAOXO  5VN1N  SVHlX  SVM  fVHXO 

0.70000  0.09799  242.0**14  0.710491751.057(7  2.00000  2.00000 
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